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Abstract 

The resummed differential thrust rate in e^e~ annihilation is calculated using Soft- 
Collinear Effective Theory (SCET). The resulting distribution in the two-jet region (T ~ 
1) is found to agree with the corresponding expression derived by the standard approach. 
A matching procedure to account for finite corrections at T < 1 is then described. There 
are two important advantages of the SCET approach. First, SCET manifests a dynamical 
seesaw scale q = i? jQ in addition to the center-of-mass energy Q and the jet mass scale 
p ~ Q\J (1 — T). Thus, the resummation of logs of "pjq can be cleanly distinguished 
from the resummation of logs of Qjp- Second, finite parts of loop amplitudes appear in 
specific places in the perturbative distribution: in the matching to the hard function, at 
the scale Q, in matching to the jet function, at the scale p, and in matching to the soft 
function, at the scale q. This allows for a consistent merger of fixed order corrections 
and resummation. In particular, the total NLO e^e~ cross section is reproduced from 
these finite parts without having to perform additional infrared regulation. 



1 Introduction 

Quantum Chromodynamics is a perturbative field theory for < 1, corresponding to energies 
above Aqcd ~ 500 MeV. However, setting up a good perturbation expansion is more difficult 
than simply working order by order in ttg. The difficulty is that when computing a quantity 
with more than one scale, logarithms of ratios of those scales appear which invalidate the naive 
perturbation expansion. For example, in e~^e~ collisions at center-of-mass energy Q, we might 
look for the distribution of jets as a function of the invariant mass p"^ of the jet. Then the 
differential cross section will have a correction of the form log^ Even for < 1 these 
large logarithms can dominate if p"^ is sufficiently small. Luckily, these logarithms appear only 
in certain combinations even at higher order, so that all the terms of the form (a^ log^ q2 )" 
be (re) summed at once. However, understanding which terms will appear, how to resum them, 
and how to combine the resummed result with fixed order results can be quite complicated. 
It is the goal of this paper to show how it can be done using effective field theory techniques. 

In this paper, resummation and matching of event shapes are studied using Soft-Collinear 
Effective Theory (SCET) [H El El IH [5] . Event shapes are observables which are sensitive to the 
overall distribution of final state particles, and therefore involve both short and long distance 
physics. We will consider mainly the event shape r = 1 — T where the thrust T is defined by 

I = max ' — — 

summing over all momentum 3- vectors pj in the event and maximizing over unit 3- vectors n. 
In the threshold region near r = large logarithms of the form log^ r appear at fixed order 
in perturbation theory. The resummed result will be valid even if as log r is large, as long as 
as and as log r are small. 

Effective field theories provide a systematic approach to resummation. They separate out 
physics at a hard underlying scale Q from physics associated with a scale of interest p^ ~ Q^t 
and from even lower scales. At each scale a separate matching calculation is done which is 
independent of physics at asymptotically lower or higher energy. Then the scale dependence 
is calculated using the renormalization group. In this way, large logarithms cannot appear 
because no two largely separated scales are accessible to the theory at the same time. The 
Soft-Collinear Effective Theory works by separating the degrees of freedom of QCD into soft 
modes and coUinear modes in different directions. The relevant scales are then associated with 
the observable of interest, such as Q^r, or with the large components of the coUinear modes. 

The resummation of event shapes in e~^e~ is not currently of extreme importance phe- 
nomenologically. It nevertheless provides a clean arena (as compared to hadron collisions) 
to explore resummation and matching. The SCET techniques and most of the formulas we 
discuss here were originally developed for 5-physics, such as resummation in 6 — > S7 de- 
cays [H El [71 [8]. They have also been applied to the study of deep-inelastic-scattering near 
X = 1 [Hllinilll] and to the production of massive jets initiated by top quark decays [I2]. One 
convenient feature of effective field theories is factorization, which allows us to use objects, 
such as soft and jet functions, calculated in one process to study another. Thus, most of the 
hard work required for calculation of the distributions we describe here can be extracted from 
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the literature. Nevertheless, there are certain features, in particular the NLO matching step, 
for which event shapes are uniquely illuminating. 

The breakdown of naive perturbation theory due to the appearance of large logarithms is 
independent of blowing up and of non-perturbative effects. To emphasize this point, and 
to simplify the resummed expressions, the running of as will be simply turned off (setting 
/3 = 0) for most of this paper. Working in the conformal limit removes one scale (Aqcd) from 
the problem and thus clarifies which large logs are being resummed. It is not hard to turn 
(3 back on, as will be shown in Section O Also, we will be including all one-loop results but 
no two-loop results. Thus, our expressions will not contain a complete next-to-leading log 
resummation, which should also resum the two-loop double logs. 

There are three steps involved in the calculation of event shapes in SCET, associated with 
three separate scales. At the hard scale ~ Q, QCD is matched onto SCET by demanding 
parton level matrix elements in the two theories be the same. This can be done at leading 
order by matching onto an operator O2 with two collinear fields, or at next-to-leading order 
by matching in addition to an operator 03 with collinear fields in three directions [131 [E] • 
The threshold resummation only involves O2 as the matrix elements of vanish near r = 0. 
(An alternative approach to C3 matching is described in At the scale /ij ~ Q\/t, the 

collinear fields freeze and can be removed from the theory by integrating them out. This 
results in the a jet function J{p^). Finally, even though we are interested in a distribution at 
the scale Qy/r, the soft degrees of freedom remain relevant down to a scale fis = Qt, after 
which they too can be integrated out of the theory. The fixed order result will have large 
logarithms of fJ^h/f^j and fJ^j/fJ^s, but in the resummed result all these logs are exponentiated 
and innocuous. 

Before we present the factorization formula and calculate the thrust distribution in SCET, 
we will review the way resummation of thrust is traditionally handled. There are a number of 
ways to resum event shapes [16l [T71 [HI [191 [201 [21], [22]. Since the focus of this work is on com- 
parisons to SCET, we will confine our attention to the original approach of Catani, Trentadue, 
Turnock and Webber [16], which will be referred to as CTTW throughout. Moreover, most of 
the other approaches reduce to [16] at next-to-leading fixed order in as (NLO) and to leading 
log, so no significant loss of generality is sustained. We will find that in the two-jet limit SCET 
also agrees with CTTW to the order we are working, although the resummed expressions are 
not exactly the same. 

A more significant difference is in the NLO matching. A critical advantage of the effective 
theory approach is that the finite parts of loop amplitudes are automatically incorporated into 
the perturbative expressions. For example, the total NLO cross section for e~^e~ is reproduced 
by combining the finite parts of the hard, jet, and soft functions and a contribution from 
a finite integral over higher-order operators in SCET. This does not necessarily entail less 
work than in calculating the total cross section through the traditional combination of real 
and virtual contributions. However, due to factorization, the infrared divergent contributions 
which are absorbed into the jet and soft functions are universal and thus the could potentially 
be used for many processes. 
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2 Perturbative QCD 



In this section, we will review some basic results from QCD on thrust, and the resummed 
expressions presented in |16j . 

To begin, consider the parton model description of e^e~ annihilation. At order the 
only process which contributes is e^e~ qq. These two jets have no structure and hence the 
cross section is simply da/dr = cro5(r), where (Tq is leading order total e+e" annihilation cross 
section. 

At order a], there are two e'^e~ — > qqg diagrams which contribute 



C^CTparton ~ 




+ 




(2) 



where the photon line on the bottom is the e"'"7'^e current coming in. The differential cross 
section is 



1 dV 

(To ds dt 



S{s)6{t) + a 



parton 



+ r + 2u 
7t 



where we have defined 



a = 



2a, 



(3) 



(4) 



and the reduced Mandlestam variables are s = {pg ^Vqf'IQ'^ and t = {pg +PqY/Q'^ with 
s + t + u = l. 

Now, for 3-parton events the thrust variable r = 1 — T reduces to 



So that 



1 da' 
ctq dr 



parton 



6{t) + a 



6{t) + a 



T = min(s, t, u). 



2(3r2 - 3r + 2) 1 - 2r _ 3(1 - 3r)(l + r) 
r(l — r) ^ r r 

-41ogT-3 



r 



(5) 

(6) 
(7) 



The second line is written to manifest the singularity structure. The remainder (ifin(T) is finite 
as r ^ 0. 

Instead of the differential thrust distribution, it is useful to work directly with the inte- 
grated quantity 



For small r, 



R 



parton 



V) ~ —2a log^ r — 3a log r. 



(9) 



Here we see explicitly the large logarithms a, log^ r and a, log r which demand resummation. 
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A useful theoretical trick, used in [TB] and in [12], is to employ a hemisphere mass defi- 
nition. This greatly simplifies the factorization formula in the two-jet limit. The hemisphere 
momentum Pl{pr) is defined to be the sum of the 4- momenta of all the particles in the hemi- 
sphere of the left (right) jet. Then thrust reduces to r = (p| -|- p|.)/Q^, and we can calculate 
it from 

da f dV pI + pI , ...s 

which we will see has a closed form expression. 

The traditional approach to resummation which we review here is due to Catani et al. [T6] . 
The basic idea is that terms of the form (alog^r)" come from multiple real collinear and 
soft emissions. The kinematics of multiple collinear emissions can be modeled by parton 
branchings, with branching probabilities proportional to the Altarelli-Parisi splitting functions 



dV 



dt dz 



da"^ X Pgg = da"- x a 



11 + z' 



t 1 



'in 



In the two-jet limit, the corrections due to soft effects can be modeled by imposing angular 
ordering. All of this follows from the coherent-branching algorithm developed in [23]. In fact, 
more general methods are available for handling soft emissions (for example, [21] )• But, to 
NLO for event shapes like thrust, the methods set up in [16] are sufficient to compare to 
SCET. 

In the two-jet limit, the factorized expression for the differential cross section in CTTW is 



dV2 



aodpidp 



R 



Jc{Q\pI)Jc{Q\pI). 



(12) 



CTTW 



The subscript on cr2 refers to the 2-jet contribution. Jq has a simple physical interpretation: 
it is the probability of finding a final state jet with invariant mass p^. Thus, it satisfies 



dp'Jc{Q\p') = l. 



(13) 



Q is a scale associated with the hard process. We will see later that the equivalent jet function 
in SCET has a different normalization. 

At leading order the jet is a massless parton and so Jc{Q^,P^) = ^{p^)- At next-to-leading 
order, the coherent branching algorithm allows one-angular ordered emission. So 



Jc{Q\p') = Kp' 



a 







q2 



1-z 



■6{p^ - z{l - z)f) + ■ 



(14) 



Here z is the energy fraction in the emission and q 



r,2 



P 



z{l-z) 
3lii 



E'^{1 — cos^). The angular 



ordering constraint for the coherence of soft emissions is implicit in the restriction q < Q; for 
> it cuts off the infrared divergences in the integral as 2; — 1. Evaluating 0141) gives 



MQ',p')=S{p') + a 



-21ogi 

p2 



+ 



(15) 
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Thus, to first order in 



1 d(j2 
(Jo dr 



5{t) + a 



J CTTW 



-4 logr - 3 



T 



(16) 



which reproduces the divergent as r ^ part of the the parton model result ([7j). 

Here we have introduced the * distribution (alternatively called the R- or //-distribution) 
which is a generalization of the +- distribution. It is uniquely defined by the two conditions 



[f{x)fr^ = fix) for a; > 



;i7) 



dx[f{x)Yr^g{x) = / dxfix) [g{x) - g{0)] (18) 
'0 Jo 

For clarity, we have added to the notation an explicit instance of the dependent variable x in 
[/ji^'^^l A useful relation is 



[/(x)]i-''^] = [/(x)]i-'n5(x) / dx'fix'). 



(19) 



This will be used extensively in the next section. 

The leading order resummation is preformed by iterating the angular-ordered emissions. 
This leads to an integro-differential equation for Jc{Q'^,p'^), similar to the evolution equation 
for parton-distribution functions. The equation is solved in an integrated form, and the 
resummed expression for the integrated thrust in the two jet limit is 



(Jo 



J CTTW 



where 



exp [—2a log^ t — 3a log r] 



r] = —2a log r. 



-2'y En 



r[2r/+ 1]' 



(20) 



(21) 

dg-2 



This equation is a combination of expressions in [T^ taken with /3 = 0. Expanding R2{t) = 
to order as the fixed order result (fT6l) is reproduced. This expression is resummed in the sense 
that R2{t) — > as r — > 0, in contrast to the fixed order result (fT6!l for da or the parton model 
expression which diverge as r — 0. 



3 Soft-Collinear Effective Theory 

Having reviewed the way the resummed thrust distribution is traditionally calculated, we 
now turn to the equivalent calculation in SCET. We will see that there are a number of 
advantages of this effective field theory treatment. Instead of a classical treatment, where 
multiple real emissions are summed at the level of the cross section, the entire resummation is 
done in SCET through the renormalization group. This makes explicit the various scales in the 
problem, and allows greater freedom to choose those scales to minimize the large logarithms 
in the distributions of interest. 
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The idea behind SCET is to separate out the quarks and gluons of QCD into soft and 
coUinear degrees of freedom. A colhnear field is associated with a hght-hke direction n^. The 
component of its momentum in that direction, p~ = n ■ p, must be much larger than any of 
the other momentum components. All the QCD degrees of freedom which could change this 
momentum have been integrated out, so p~ appears as a fixed label. For example, a colhnear 
quark is written as C,n{p) with a fixed p~ . In addition to colhnear quarks and gluons, there 
are soft quarks and gluons. These soft fields can only interact with each other or transfer 
momentum to the soft components p~^ = n ■ p oi a colhnear field. 

The Lagrangian of SCET has a separate gauge invariance associated with soft gluons and 
gluons in each colhnear direction. It is useful to maintain this gauge invariance explicitly in 
the operators of the theory with the use of Wilson lines. For example, a two-jet operator is 



where 



O2 = ^nWnYnYY^WlU, (22) 

Wn{x) = PexpS^ig j dsn ■ Anifis + x)| (23) 

Yn{x) = Pexp i^ig J dsn ■ As{ns + x)| . (24) 

The colhnear Wilson lines Wn maintain colhnear gauge invariance and the soft Wilson lines 
Yn maintain soft gauge invariance. 

The starting point in the effective theory approach to event shapes is again factorization 
in the two-jet limit. In SCET, the event-shape distributions near the endpoint come from 
matrix elements of O2 [121 [3 [23 [2S] • In terms of the hemisphere masses pi and pn defined 
above, factorization implies 



dV2 



dpi dp% 



SCET 



|C^^(/i) p J dkLdknJipl - QkL, ^i)J{pI - Qkn, fi)S{kL, kn, /i). (25) 



The hard function Ch{^)i the jet functions J{p'^,fi) and the soft function S{k,fi) will be 
defined below. 

The form of this factorized expression fl25|) has a physical explanation. Each jet function 
J(p^,yu) comes from one of the colhnear quarks. It represents, like Jc, the probability for 
producing a jet of invariant mass from that colhnear field (the precise relation to Jq is 
explored below). Recall that the large component of the momentum of a colhnear field p~ = 
n-p is fixed. Since, in the two-jet limit, the jets are back to back with center of mass energy Q, 
we must have p~ ^ Q- The small component p'^ = n ■ p can vary. If there are no soft effects, 
then the hemisphere mass is simply p'jj = p^ = p~p'^. However, as the factorization formula 
implies, the colhnear jet can give up some soft momentum to the soft QCD background, 
leading to p'^ ^ p^ — k. The hemisphere mass is unchanged by this emission, but now 
p^ = p~{j)^ — k) = pIj — Qk, which explains the form of fl25l) . 

The function Ch is a hard function, which comes from integrating out hard modes of QCD 
in matching to SCET. Demanding 

{qi'q) = Ch i^n^^n) (26) 
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for states with two external quarks lets us calculate Ch order by order in perturbation theory. 
The computation entails taking the difference between QCD and SCET graphs, such as 



C 



H 





(27) 



where the ® refers to an insertion of O2, and the right diagram is only a representative 
contribution (for details, see [I3]). The difference is finite because the infrared divergences in 
QCD and SCET are the same and the UV divergences are removed with counterterms. At 
one-loop the matching gives [9l [13] 



IH log 



with Ch 



a 



Th = —ot and 



The anomalous dimension and renormalization group evolution of Cni^J') can be found 
from the UV divergences of the same 1-loop graphs. The RG equation is 



with solution 



dlog fi 



C/^(/i) = C^(/i,j)exp 



-2r^,iog\ 



27/f)CH(/i), 



(29) 



. 2 /i^ 



1 

IH log — 



-r„iogiiy 



(30) 



To first order in the dependence drops out and fl28l) is reproduced. The natural matching 
scale /i/j = Q is manifest in this expression. 

The jet function J(p^,/i) is the imaginary part of the propagator of a coUinear quark. It 
is defined by [271 [12] 



J(p^/i) 



Tiin ■ p) 



Im 



(31) 



It can be thought of as the spectral density for a jet of coUimated particles interacting with 
a soft QCD background. The jet function can be calculated order-by-order in perturbation 
theory through the discontinuity of conventional Feynman diagrams, such as 



J(p^, /i) ~ Disc 



+ 



+ 



+ 



(32) 

Some of the infrared divergences in these graphs are cut off because the coUinear quark is taken 
to have invariant mass . These calculations have been done in the context of 6 57 [6], [7] 
and for deep-inelastic scattering [Qj [TT] and the jet function is known to two loops [27]. Due 
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to factorization, the same jet function applies in these processes and in e~^e annihilation. To 



first order in a^, the result is 



J(p^/x) = 5(/)[l+c,] + 



rjiogg, + 7j 

p2 



(33) 



with cj = a ^1 — ^ j , Tj = 2a and 7j = — fa. 

The renormalization group evolution of the jet function, in contrast to that of the hard 
function, is non-local in p"^ 



dJ(p^/i) 

dlog/i 



-2rjlog^-27 



J(p,/x) + 2rj/ . 



(34) 



This is similar to the Altarelli-Parisi equation for the evolution of the parton distribution 
functions (pdfs). In contrast to pdfs, however, the jet function J(p^,yu) is perturbative as long 
as > Aqcd- Simplification is achieved through use of the Laplace transform 



(35) 



whereby the evolution becomes local in u 



-2Tjlog- 



1 



d log yU \^ ^ n'^ue'^^ 
Now the RGE can be solved like that of Ch'- 



27j 



(36) 



]{v,^^L) = exp 



with 



— log — -7jlog — 
. 2 /i,' f^l 



Vj = -rjiog-2. 



(37) 



(3^ 



Finally, the inverse Laplace transform produces [TT] 



J(/,/i) = exp 



— log — - 7j log — 
. 2 /^l 



2 \ 



where 



(39) 



(40) 



The functional dependence on dr^^ in j{drj.) comes from a functional dependence on log ^ in 
the fixed order expression for the jet function (133|) at the matching scale fi = /ij. 



8 



Finally, we have to calculate the soft function S^k^, kji, It is defined though matrix 
elements of soft Wilson lines [251 W2\ 



s{kL,kR,fi) = Y,m^yn\x){x\YM*\o). 



(41) 



X 



Unlike the jet function, the soft function is often evaluated at scales k ~ Aqcd, where it 
is non-perturbative. It can, in principle, be measured experimentally and then be evolved 
perturbatively, like the pdfs. Or it can be modeled [221 EOl [12] • However, in our case, since 
we are taking to be small and fixed, it can be calculated in perturbation theory. Even if S 
is non-perturbative, the perturbative calculation is useful because the UV divergences dictate 
the anomalous dimensions and hence the evolution equation. 

The soft function, being completely defined in terms of Wilson lines, can be calculated in 
QCD through diagrams such as 



S{kL, kR,fi) 




+ 




+ 




+ ■ ■ 



(42) 



where the kinked lines are Wilson lines in the n and n directions. To order as, S{kL, k^) can 
be written as a product 

S{kL,kR,n) = S{kL,fi)S{kR,fi) (43) 



where [281 



S{k,ix)=5ik) [l + cs] + 



Ts log ^ + 75 
k 



(44) 



2 

with cs = «Y2! = —4a and 75 = 0. 

The evolution of the soft function involves the same kind of non-local equation as for the 
jet function. This can be seen from examining its anomalous dimension directly, or from RG 
invariance of the convolution appearing in the factorization formula fl25|) . The solution is 
therefore similar. Explicitly, 



S{k, fj,) = exp 



— log 7s log — 

2 /is /is 



s{dn 



1 ''^ 
k \fis 



where 



and 



TT 



2 r 

12 2 ''^ 



Vs 



-Tq loe 



/is 



(45) 

(46) 
(47) 



Expanding to first order in ag reproduces ( 144|) . 
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Now we have all the ingredients appearing in the the factorization formula fl25l) . At order 



a,, the result is 



1 dVa 
(Todpidpl 



+ a6{pj 



2 log 



SCET 

3' 

2 



1 + a 



pI 



[pi 



+ a5{p\ 



2 log 



P%. 



TX 



(4J 



Note that the dependence on all the scales /i, /i^, /ij and /i^ has completely canceled. This is 
a non-trivial result which requires three relations among the six anomalous dimensions: 



iH + lJ + 



^ = 
2 



(49) 



and is a strong consistency check on the entire formalism. 

Already SCET can be compared to CTTW. The differential distribution in SCET, Eq. ( l48i) . 
matches the CTTW distribution, Eq. (fT2|) . when the jet functions are expanded to first order 
using Eq. (fTSjl . The SCET expression has an additional finite piece (the — 1 + ^ term) which 
comes from loop graphs which do not enter into the traditional formulation. 

Still working at leading order in a^, SCET produces a simple form for the thrust distribu- 
tion near r = 0: 



1 dcr2 
(To dr 



5{t) 



J SCET 



a 



a 



-4 logr - 3 



(50) 



This reproduces the leading behavior for small r of the both the parton model expression, 
Eq. (j7]), and the resummed expression with the functions Jq expanded to first order, Eq. (ITBi) . 

To get the resummed thrust distribution from SCET, we need to calculate a couple of 
convolution integrals. First, the soft and jet functions must be combined into 



K{p'^,n) = J dkJ^p"^ — kQ, jj)S{k, n). 



(51) 



Second, the K functions must be integrated against the event shape. For thrust, we need 

R2{t) = \Cn{l^)? j Kipl^^)K{p%^^)eir-^4A)■ (52) 



Both of these convolutions can be evaluated by performing the Laplace transform. 

Note that the function names can be misleading. Here, Ch{p-)K{p'^, fi) (and not just the 
SCET jet function J{p'^,fi)) is playing the role of Jc(p^) from the CTTW formulation. Thus, 
the CTTW jet function is a combination of the soft and jet functions in SCET, as expected 
because the coherent branching algorithm used to derive it incorporates both soft and collinear 
effects. 



^ The relations among the Fs can be understood on more general grounds by relating the Fs to a universal 
cusp anomalous dimension |31| . 
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Solving the K function with the Laplace transform techniques gives 

,2 ,,21 rr _ ,, ,, n /g^ 

■ ^iug 75iug — I 



K{p^,fi) = exp 



— log — - 7j log — 

2 ;,2J 



where 



and 



X 



— log - log —^dr,, + k{d^. 



2 \ 



1 + Cj + Cs 



TT 

12 



(rj + r< 



J LP V/^j 



dl + {ij + is) d, 



Vk = Vs + Vj = 2a log— . 



(53) 
(54) 
(55) 
(56) 



Note the residual logarithms of Qfi^/fi'^ in line f l5^ which are not resummed. These imply 
that we cannot choose /ij = fig- Instead, the natural scales (those which remove the large logs) 
should satisfy fig = fi'j/Q- Then, we can evolve either the jet function from fij down to fis or 
the soft function from fig up to fij but there is not a single scale which can minimize all the 
large logs. 

Finally, to evaluate the convolution for the thrust integral, we use the Laplace transform 
of the ^-function 

1 



up 



(57) 



Combining this with the expression for Klp"^, fi) above, the SCET prediction for thrust in the 
2- jet limit is 



J SCET 



2alog2 4 + 3alog4-«log'4 
/i| /ij /il 



/i^ 



-4o lO! 



X 



-2a log^ + 4a log ^^d2n^ + k{d2r^^, 



Q}is 



2 /o2^\ 2r?fc -27^% 



r[2r/fc + 1] 



(5^ 



Note that this expression is explicitly independent of /i. 

Now, choosing the natural scales /i^ = Q, \ij = Qy/r, and Hs = ji^Q = Qt to remove the 
large logs, the thrust distribution becomes 



-R2{r) 
ctq 



J SCET 



exp [—2a log^ r — 3a log r] r (9^) 



r[2r/ + l] 



with 



TT _ 

— la a 

3 ' 3 



-a9^ — -a(9„ 
2 2 ^ 



(59) 



(60) 



and 7] = —2a log r as in Eq. (12T1) . 

This expression would be identical to the CTTW expression, Eq. (120|) . if f = 1. Recall that 
the function f comes from 1-loop matching in SCET, which turns into boundary conditions 
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for the renormalization group evolution. In the approach of CTTW the boundary condition 
is simply that Jc{p^) = ^{p^)- Nevertheless, the effect of f 7^ 1 to order is only to provide 
a finite constant, as can been seen from comparing equations fl50p and fll6p . That is, there 
is no contribution of order rj in the difference, only of order 77^ and higher. Beyond order 
as, the dri does change the rj dependence. But since r] ~ a^logr these terms are subleading 
to the dominant log^ r and can get corrections from higher-loop effects. For example, at 
two-loops, a term of the form log^ r is generated. Thus, SCET and CTTW agree to first 
order in and for the resummation of the leading large logarithms, which is the order to we 
have been working. 



4 Matching to hard emissions 

Now let us look at how fixed order results and resummation are combined. There are two 
elements to this: (1) matching to hard emissions to get the differential distribution correct 
away from the two-jet region (i.e. away from r = 0); and (2) including finite parts of loops to 
reproduce fixed order inclusive results. 

Let us begin with the matching procedure described in To match to the hard emissions 
at order Og we need the parton model differential cross section from equation ([7]). The divergent 
part is already contained in the two jet contribution, so the remainder is 

DUr) = [ dr'dUr') (61) 
Jo 

= logr[6r + 41og(l - r)] - 21og'(l - r) + 3(1 - 2r) log(l - 2r) - 4Li2(-^) + 6r + ^r^. 

1 — r 2 

The final result for the matched integrated thrust distribution from [16] is 



—R{t) = (1 + (7i) <^ (1 + Ci) exp [-2a log' r - 3a log r] + aD^r^ir) \ , 

.^0 J CTTW L i [^?7 + ij J 

. (62) 

where ai = |a = ^ is the NLO contribution to the total cross section and Ci = a — |j 

is chosen so that R{\) = 1 + ai (recall that r < at order a^). In the QCD approach 
the factor cxi = — must be be determined by a independent calculation of the total rate. 
This involves combining the real and virtual contributions at order using an appropriate 
infrared regulator. As emphasized in [TB] there is arbitrariness in the matching because the 
hard emissions (in -Dfin(T)) are fixed order but the two-jet contribution is resummed. For 
example, to the same accuracy -D(r) could be multiplied by the exponential. 

Now, let us turn to the matching in SCET. To include thrust distributions away from the 
endpoint, we can either attempt to add power corrections to SCET, or we can match to higher 
order operators as described in [131 [S]. Matching is much simpler. To perform the matching, 
we add new operators 
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where the • • • are the additional terms coming from Wilson lines necessary for gauge invariance. 
A!^^ is a coUinear gluon in direction 77,2- The matching demands that 



+ = (QCD) 



(64) 



where the subscript means the matching is done at the hard scale fih- 

There is some arbitrariness in the definition of the matrix elements in SCET due to repa- 
rameterization invariance. The matrix element of a coUinear quark C,n{p) on a QCD quark 
state \q) is only defined up to its soft momentum component n ■ p. Moreover, a basis for 
summing over directions n for O2 and Ui, n2, and for must be chosen as well. A certain 
convention was described in [13] for resolving these ambiguities and others are possible. In 
any case, while different conventions may shift the contributions from O2 and in fl6^ . the 
sum is parameterization invariant. Thus, independently of the convention we have 



[dcrl'^'^ 



J SCET 



[da] 



parton ; 



(65) 



where [dajparton ~ (QCD)^ is the tree-level parton model cross section, as shown in Eq. 

Now, we already know that {02)^ gives fl50l) to first order in and that the parton model 
distribution at order is ([7]). Thus, with obvious implicit phase space factors. 



{02 + Os)i - ml = d 



filT 



(66) 



So, at leading order, the contribution from the sum of (Oa)^ and the interference between 
O2 and to the differential cross section is unambiguous. The running of could also be 
included even though it does not resum any large logs for the event shapes under consideration. 
With running, the matrix elements of O2 and would appear with separate renormalization 
kernels, and so the final differential cross section would end up depending on the conventions 
chosen. The ambiguity could be resolved by a careful higher order treatment, but for the 
purposes of comparing to the CTTW prediction for thrust, we will simply not renormalize the 
finite terms. 
Then, 



-R{t) 



Q-2-fEri 

= exp [-2alogV - 3a log r] f{dn) 

SCET ■'-[-'- + -^Vl 



(67) 



With r{drj) and 77 as in Eq. (160|) and Dh^It) in fl6T]) . The total cross section in SCET is given 



SCET 



1 + a 



TX 



aD 



fin I 



TT 



(6^ 



This is the correct total e"''e~ total cross section to first order in a^! 

Let us review the contributions that go into the cross section. First 
Hh = Q, there is the finite part of the loop matching to QCD, jc/fp = 1 + a 
1 + 2.3cri. Next, at the scale fi"^ = = Q'^t where we integrate out the coUinear fields, the 



at the hard scale 



jet functions give 2cj = a{7 



TT 



-l.QcTi. At the seesaw scale /i^ = p'^/Q, the soft function 



gives 2cs = IT = l-lo"i, and finally the finite part of the real emission integral, away from 
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0.0 0.1 0.2 0.3 0.4 0.5 

1-T 

Figure 1: Matched resummed differential thrust distribution in SCET and in the standard 
approach of CTTW with fixed couphng as = 0.4 



r = gives 0"3 = a ^| — = — O.Scri. In producing the total cross section, only the soft 
and jet contributions are infrared divergent. However, their convolution, which appears in the 
function K{p^,^) is infrared finite. The hard matching and the hard emissions are IR finite 
by themselves. Thus, the total cross section can be understood as a combination of a process 
dependent IR finite hard part and universal but IR regulator dependent contributions from 
soft and collinear emissions. 

The SCET thrust distribution ( 167|) is compared to the CTTW thrust distribution ( 162|) 
in Figure [H The plot is of ^ = R'{t) with a, fixed at 0.4. A more careful rendition of the 
differential thrust distribution would take the derivative of R{t) before assigning the matching 
scale (ij = Q^pT^ as is done below. However, the effect is higher order, and so we just plot 
i?'(r) directly. 

5 Generalizations 

In this section, some simple generalizations are described. The above results were derived 
assuming to be constant in order to emphasize the resummation of Sudakov logarithms in 
contrast to large logarithms associated with the scale Aqcd- Now, it is shown how the results 
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can be modified with running a^. Also, the SCET prediction for another event shape, the jet 
mass p, is given. 

It is straightforward to allow to run. Including 1-loop running, the effect is to modify 
the single and double logs in the following way |llj : 



-1 2 

a log — 

V 



a log 



S'(z/,/i) 



47rCf 



log 



For example, the differential thrust distribution for r > with running becomes 



1 da 
(To dr 



J SCET 



^ exp [45(g, gr) + 6A(g, gr) - 8SiQ^/T, Qr) - 6A{Qy/T, Qt)] 
+— 4n(r), 



where 

dfin is given in Eq. ([7]) and 

2 

+ 



r/ = 4A(gv/7, gr) 



^)«,(g) + (7 - ^)as{QV^) - y«,(gr) 



(69) 
(70) 

(71) 
(72) 



'-as{Q\/T) - as{Qi 



(73) 



This is the same function f as in Eq. (!60l) . but with the factors evaluated at the appropriate 
matching scales. 

At this point, one would hope to compare to data. However, besides the obvious shortcom- 
ing of not containing the full NLL resummation (it does not include effects of the two-loop 
cusp anomalous dimension), this parton- level expression does not include important non- 
perturbative effects. Due to the running of as, the perturbative expression breaks down when 
the soft scale is of order Aqcd, that is, when r ~ Aqcb/Q, as can be seen exphcitly in fl73l) . 
In fact, even for significantly larger values of thrust power corrections of order Aqqb/Q be- 
come quantitatively important, at least at LEP energies. This problem has been approached 
elsewhere using SCET [251 ES] and with other techniques [211 [23 ED]- 

Other event shapes can be studied the same way as thrust. For example, consider the 
heavy jet mass p defined by 

(74) 



In this case, the matching scales are ph = Q, Pj = Qy/p and ps = Qp and SCET gives for 
p > 



1 da' 
ao dp 



J SCET 



P 
X 



ch\^ exp [45(g, Qp) + 6A(g, Qp) - %S{Q^p, Qp) - 6A{Q^, Qp)] 



k{dr, 



~1eV 



~k{d,) 



2a, 



d&n{p), 



(75) 
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where 

r, = 4A{Q^,Q). (76) 

This formula agrees with the jet mass distribution in [16] to leading log and first order in a^. 
The same function appears for jet mass and for thrust because to order in the parton 
model, p = T = min(s, t, u). 

6 Conclusions 

We have investigated how to combine resummation with next-to-leading order matching of 
event shapes in the original approach of [16] (CTTW) and using SCET. In the CTTW for- 
mulation, real emissions from coUinear splitting functions are used and various kinematical 
features associated with soft emission are combined to derive a differential cross section. The 
cross section factorizes into the product of two jet functions. Resummation is done by solving 
a differential equation for the jet functions in terms of the physical scales and of the 
event. In contrast, SCET factorizes the event shape distribution into a contribution from 
hard, jet, and soft functions. These functions are matched at the scales ph = Qi f^j = P and 

2 

/is = ^ respectively and resummation is done through renormalization group evolution. The 
resummed thrust distribution in SCET and CTTW have been compared, and found to agree 
to next-to-leading order in and for leading-log resummation. 

The resummation of thrust in SCET brings to light a number of features of the process 
not obvious in CTTW formulation. For example, the appearance of the seesaw scale fig = 
Q{1 — T) as the natural matching scale for the soft function follows from the kinematics of 
the SCET factorization theorem. Of course, the existence of this scale has been known for 
a long time from QCD, but in the effective field theory approach this scale just drops out of 
the factorized expression. Thus, for more complicated processes, it is reasonable to expect a 
similar transparency for the matching scales, which may facilitate subleading log resummation. 
In fact, two- loop, three- loop, and some four- loop anomalous dimensions for various quantities 
are already available [32], [331 IMl [27] , and so subleading log resummation appears quite feasible. 
The biggest impediment to using these more accurate resummed results in a comparison to 
data is that power corrections of order Aqcd/Q have an important quantitative effect on event 
shapes. However, these corrections should modify only the soft function while higher order 
resummed expressions for the the hard and jet functions will remain valid. Thus, the effective 
theory should be able to weave together the perturbative and non-perturbative calculations. 

A new result of this paper is the demonstration that inclusive quantities, such as the total 
cross section for e~^e~ can be calculated in a new way using SCET. Instead of taking the full 
differential n+1 body cross section and combining with the one-loop n-body cross section, 
SCET combines finite parts of loops of soft and jet functions with a hard matching calculation 
and a finite integral over hard emissions. The soft and jet functions depend on the infrared 
regulator, but their convolution, and the hard function, do not. For e+e" annihilation at 
NLO, this may not be so impressive, but the procedure promises to apply to more complicated 
processes, perhaps even some for which NLO results are not available. It would also be very 
interesting to explore this mechanism at NNLO or to work with hadronic processes where the 
singularity structure is more complicated. 
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